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Abstract 



We analyze possible implementations of quantum algorithms in a system 
of (macroscopic) Josephson charge qubits. System layout and parameters to 
realize the Deutsch algorithm with up to three qubits are provided. Special 
attention is paid to the necessity of entangled states in the various implemen- 
tations. Further, we demonstrate explicitely that the gates to implement the 
Bernstein- Vazirani algorithm can be realized by using a system of uncoupled 
qubits. 
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Quantum information processing has initiated an impressive research activity throughout 
the scientific community during the past decade. The interest in quantum computation, in 
particular, is stimulated by the discovery of quantum algorithms 0,0] which can outperform 
their classical counterparts in solving problems of significant practical relevance. 

In recent years, considerable progress has been made in the field of "quantum hardware" , 
i.e. in the search for physical systems appropriate for the implementation of quantum com- 
putation. It is equally important in the developement of quantum computing to experimen- 
tally realize complete quantum algorithms. To date, this has been achieved in liquid-state 
NMR l^ini, in atomic physics |]TT|] and by optical interferometry |]T^. A solid-state imple- 
mentation of Grover's algorithm has recently been proposed in Ref. 
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The quest for large scale integrability has stimulated an increasing interest in supercon- 
ducting nanocircuits [|I4|-|T8| as possible candidates for the implementation of a quantum 
computer. The recent experimental breakthrough for Josephson qubits [p!9|-|2l[| is the first 
important step towards a solid-state quantum computer. 

Naturally the question arises whether, at the present technological level, it is possible 
to implement also quantum algorithms in these systems. Here we concentrate on charge 
qubits ||14|-0 and show how the Deutsch algorithm p^-|25| and the Bernstein- Vazirani al- 
gorithm |]2^ can be run on a Josephson quantum computer. We analyze the experiment by 
Nakamura et al. in terms of quantum interferometry |2|] and show that it corresponds to 
the implementation of the one-qubit version of Deutsch's algorithm. By generalizing this idea 
we show how the A^-qubit Deutsch algorithm, with < 3, can be implemented. Entangle- 
ment is required only for A^ > 3 . Finally we show explicitely that the Bernstein- Vazirani 
algorithm |^ can be implemented using uncoupled qubits (for arbitrary A^). Therefore it 
can be realized by means of the setup of Ref. [|T9[ . 

Consider the subset of Boolean functions / : {0, 1}^ — > {0, 1} with the property that / is 
either constant or balanced (that is, it has an equal number of outputs as Is). The Deutsch- 
Jozsa algorithm P^|-p^ determines whether a function / is constant or balanced. With a 
classical algorithm, this problem would, in the worst case, require 2^~^ + 1 evaluations of / 
whereas the quantum algorithm solves it with a single evaluation by means of the following 
steps (here we focus on the refined version by Collins et al. |2^, see also Fig. 1). 

(i) All qubits are prepared in the initial state |0), therefore the A^-qubit register is in the 
state |00. . .0). 

(a) Perform an A^-qubit Hadamard transformation Ti 
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where x-y = (xi A © . . . © (xiv A i/n) is the scalar product modulo two. This is equivalent 
to performing a one-bit Hadamard transformation to each qubit individually. 
(Hi) Apply the /-controlled phase shift Uf 
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Note that we will use the convention /(OO ... 0) = 0. 

(iv) Perform another Hadamard transformation Ti. 

(v) Measure the final state of the register. If the result is |00. . .0) the function / is 
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constant; if, however, the amphtude a|oo...o> of the state |00. . . 0) is zero the function / is 
balanced. This is because 
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Entanglement is not needed for the one-bit and two-bit case p5| while it is necessary for the 
exponential speedup for a higher number of qubits ||26l . 



In order to implement the algorithm we have to show that each individual step (prepa- 
ration of the state, gate operations, measurement) can be realized. It is well known how 
to prepare and to measure the states in Josephson charge qubits |]14|,|16|,|19[ . Our task is to 
demonstrate that the gate operations corresponding to all possible functions / can be per- 
formed with a single device. An important aspect of our proposal is that the gate operations 
are represented in a basis of superpositions of charge states. 

One-qubit and two-qubit Deutsch algorithm - In the one-bit case there is one constant 
function and one balanced function / (due to our choice /(O) = 0, see above). The gate 
Uf implementing the constant function is the one-qubit identity operator Ii. The balanced 
function can be represented (with respect to the computational basis) by cr^ where cjj denote 
the Pauli matrices. 

The sequence of steps (i)-(v) can be carried out with a Josephson qubit. A Josephson 
charge qubit |l^Jl6[] is a Cooper-pair box (see Fig. 3a) which is characterized by two energy 
scales, the charging energy Ech = (2e)^/(2C) (here C is the total capacitance of the island) 
and the Josephson energy Ej -C E^h of the tunnel junction. At low temperatures T <C Ej 
only the two charges states with and 1 excess Cooper pair on the island are important and 
the system behaves as a two- level system with the Hilbert space {|0), |1)} and the one-qubit 
Hamiltonian 



= (Eeh/2) (2n,, - 1) a, - {Ej/2) a, . (4) 

Here rix = CgVg/{2e) is the offset charge which can be controlled by the gate voltage. 

The one-bit version of the Deutsch algorithm is already realized in the experiment by 
Nakamura et al (see Figs. 1 and 2). First the system is prepared in a symmetric superposition 
of the states. The Rabi oscillation in the experiment corresponds to the action of the 
controlled phase shift Uf. Finally the system is measured. Note, however, an important 
difference between the experiment and the steps (i) - (v): In the usual representation of 
Deutsch's algorithm the gate Uf acts on the same basis states which then are measured; 
the Hadamard transformation produces the appropriate superpositions. In contrast to this, 
the measurement in the experiment is done in the charge basis while the "gate" acts in 
the basis {|+), |— )} which is related to the charge basis by a Hadamard transformation. 
By identifying Uf <-> exp {i{Ejt/2h)az) the different functions / can be implemented by 
choosing the time t appropriately (see Table I). 

This observation suggests the possibihty to implement the Deutsch algorithm in a setup 
of more than one charge qubit by performing the same sequence of gate voltage pulses as in 
the experiment of Ref. [T^. The gates should operate at the degeneracy point n^^^ = 1/2 of 
the charge qubits. What we need is to find the proper parameters and operation times to 
obtain all possible gates Uf. 
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For two qubits this realization is obvious as the two-qubit algorithm can be implemented 
by using two uncoupled qubits ||25|. The gates a^}^ I'l'^ ®<yf\ u^}^ ®crp^ implementing 

the balanced functions (the upper index denotes the qubit number) and for the 

constant function can be realized in complete analogy to the one-qubit algorithm. 

Three-quhit Deutsch algorithm - The realization of the three-qubit version of the al- 
gorithm is more difficult. Apart from the constant function 35 balanced functions need 
to be implemented. Moreover, the three-qubit algorithm involves gates Uf which produce 
entangled final states. 

The goal is to proceed along the same lines as above, that is, preparation of the state 
1 000), sudden sweep of n^^^ etc. The action of the gates Uf takes place in the basis 

{| + ++),! + H — ), • • • , I )}• In order to find efficient ways for the implementation 

we first analyze the functions / and the corresponding gates Uf. 

Apart from the constant function and its gate there are 7 balanced 

functions for which the gates are separable: cr^^^ Ii'^ Ii'^ ® cr^^-* ® /j^'*, . . . , cr^^^ ® 
crp-* (J^p. Further there are 12 balanced functions for which the gates factorize into a 
one-qubit part and a two-qubit part as in example /) below. The other gates entangle all 
three qubits and can be divided into two classes (see example II) and III))- There are 12 
gates of class II) and 4 gates of class ///). 

/) i {l['^ ® + a« ® - ) ® + a« ® a^)) ® a^^^) 

//) ^ ® ) ® Jf ) - ® ) ® af^ + ® ® if + ^ ® of ® of) 

III) \ (of ® if ® if - if ® of ® if + ® if ® of + of ® of ® )) 

All separable qubit operations can be carried out in analogy with the one-qubit case above. 
In the following we discuss how the entangling gate operations can be achieved. For the 
reahzation of these gates a coupling of tunable strength between the qubits is required. 

There are various ways to couple charge qubits |l^ - pJ|jr8| . Here we investigate coupling 
via Josephson junctions |28|. Each qubit island is coupled to its nearest neighbor using a 



symmetric SQUID (see Fig. 3b). 

Assuming that both the j-th qubit and the j-th coupling junction are tunable by local 
fiuxes the Hamiltonian for the iV-qubit system at the degeneracy point nf = 1/2 

reads 



N 

r(j + l) 



i=i (5) 
-(1/2)J^')(<I>^)) + h.c. 



Here jf is the Josephson energy of the j-th coupling SQUID and o± = {o^ ± iOy)/2. 
For small coupling capacitance Cf < C^^^ we have Ef = {Cf /C^^^)Ef /2. We will 
assume that Ef is negligible. Since in practice the capacitive coupling is always present it 
is necessary to have Jk\^ = 0) ^ AEf . Then the dynamics of the system approximates 
the ideal dynamics sufficiently well. 

Consider now the first and the second qubit coupled by jf. By choosing, e.g., —Ef = 
Ef = ±jf = J and the operation time t ~ 0.97(27r/J) we obtain an operation similar 
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to a swap gate for the qubits 1 and 2 for which we introduce the notation (in the basis 

{|++),---,|--)}l2) 



[±12] :=: 



f ±i\ 
10 
10 
V±i / 



(6) 



By denoting one-bit phase shifts for the j-th qubit 

'l \ 



±t 



(7) 



we can write a sequence of operations which gives the two-bit entanghng gate in example I) 
above: 



[+l][-2]-[-12]-a 



(3) 



(8) 



After suddenly sweeping n^^ and n^^^ to the degeneracy, first the one-bit phase shifts are 
performed while J^^ = 0. Then J^^ is switched on suddenly in order to do the two-bit 
rotation. The rotation can be done at any moment since the third qubit is decoupled 
from the other two. Finally the n^/^ are swept back suddenly and the register is measured. 

There are numerous ways to represent the three-bit entangling gates. At least two 
different two-bit rotations need to be applied. During the second two-bit rotation the third 
qubit has to be "halted" . This can be done by switching off both the Ej and the Jk which 
couple to this qubit. A possible sequence for example II) is 

[+l][-2]-[+13]-[-12] , (9) 



and for example ///) 

^i'^^?Vf-[+12]-[+23]-[+12] . (10) 



The complete set of entangling gates can be obtained from the sequences (^ - (|r^) by cyclic 
permutations of qubit numbers (and appropriate sign changes), thereby paying attention 
that the parameter settings are compatible for both one-bit and two-bit operations. We note 
that universality in quantum computation implies that given some two-bit gate, it is possible 
to realize any algorithm. This procedure, however, in general requires much longer sequences 
of one-bit and two-bit gates. Our proposal allows the implementation of simple algorithms 
following the scheme given in Fig. 1 and is amenable of an experimental verification with 
present-day technology. It is interesting to note that the completely entangling gates of class 
//) and ///) can be realized approximately with a single three-qubit operation. In Table II 
we list the parameters for the various implementations including estimates for the accuracy 
of the respective operation. 

The Bernstein- Vazirani algorithm P7|J2^ - It is analogous to the Deutsch algorithm 
described in the beginning, with the difference that the function / has the form f = a- x(Bb, 
{a,x G {0,1}^, 6 G {0,1}). By measuring the register after running the algorithm once 
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(the gate in step (in) is denoted by Ua) one obtains the number a in binary representation 
which classically would require N function calls. The fact that there is no entanglement 
in the Bernstein- Vazirani algorithm has been observed in Ref. Here we demonstrate 



explicitely that the algorithm can be obtained by applying only one-qubit operations. 

The gate Ua can be rewritten as a product of single-qubit gates (using the definition 
{akT ■■=■■ h) 

N 

Ua = {-lf Ui^^^P (11) 

i=i 

where aj denotes the j-th digit of a in binary representation. Apart from the global phase 
(— 1)^ this is the part of the Deutsch algorithm with completely separable gates. As the 
algorithm starts with a product state, no entanglement is involved at any step. (We note 
that one can rewrite the action of the whole algorithm TiUa'H as Y{j{(T^''Y^ which trivially 
gives the result.) It is therefore possible to realize the Bernstein- Vazirani algorithm with 
Josephson networks in complete analogy with the implementation for the one-qubit and 
two-qubit Deutsch algorithm. 
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TABLES 
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time t 


constant 
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27:%/ E J 
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TABLE L 



gate 


implementation 


Ef 


Ef 


Ef 




7(2) 


7(3) 


time t/{2'Kh/}) 


a|o...o) 
{E^^ = 0) 


«|0...0) 

{e'^^ = J/40) 


II) 


sequence (9) 


-J 


J 


J 


-J 
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0.97 (2bit op.) 


2 • 10-3 


2 • lO-'' 


II) 


single operation 


-J/2 





J/2 


J 


J 





0.80 


7 • 10-5 


2 • 10-2 


III) 


sequence (10) 


J 


-J 


J 


J 


J 





0.97 (2bit op.) 


3 • lO-'' 


6 • 10-3 


III) 


single operation 


J/2 


-J/2 


0.83J 





J 


J 


1.19 


< 10-5 


2 • 10-3 



TABLE IL Parameters for various realizations of the gates II) and ///). The coefficient a|o...o) 
is a measure for the fidelity of the operation (for an ideal operation a.|o...o) = 0)- The operation 
time for the sequences refers to the time needed for the two-qubit rotations. Single-qubit rotations 
are assumed to be perfect. 
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FIGURES 



FIG. 1. The sequence of operations to perform the Deutsch algorithm on a register of N qubits. 
According to Ref. it can be interpreted in terms of quantum interferometry. The first Hadamard 
transformation produces a superposition of all possible states. Thus, with the application of the 
/-controlled gate Uf the outcome of / for all possible arguments is evaluated at the same time. 
The second Hadamard transformation brings all computational paths together. 



FIG. 2. Schematic representation of Nakamura's experiment |19]. The qubit is prepared in 
the ground state |0). After suddenly sweeping the gate voltage the system starts Rabi oscillations 
between the eigenstates of the new Hamiltonian |ib) = (|0) it |l))/\/2. After the time t the gate 
voltage is swept back suddenly which freezes the final state; then the qubit is measured. 



FIG. 3. a) A charge qubit. The Josephson energy of the junction can be controlled by the 
magnetic flux $: = 2£j cos (tt^/^q), where £j is the Josephson energy of the junctions of the 

symmetric SQUID and <&o = h/{2e) [|16|. Typical parameters are <fj ~ 30//eV and -Ech ~ 500^eV. 
b) A possible realization of coupled charge qubits. 
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